Possible superconducting symmetry on doped J1-J2 model 
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By making use of renormalized mean-field theory, we investigate possible superconducting sym- 
metries in the ground states of ti-tz-Ji-J-z model on square lattice. The superconducting symmetries 
of the ground states are determined by the frustration amplitude t%/t\ and doping concentration. 
The phase diagram of this system in frustration-doping plane is given. The order of the phase 
transitions among these different superconducting symmetry states of the system is discussed. 
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Frustrated magnets based on transition metals have at- 
tracted much theoretical and experimental effort in the 
past years, a variety of exotic quantum effects 1 -' 2 -' 3 -'^ ex- 
pected when competitive interactions lead a system into 
a frustrated state, where it is impossible to satisfy all the 
pair interactions simultaneously. One of the most simple 
yet very important systems is spin-1/2 Heisenberg model 
on a square lattice with competing nearest J\ and next- 
nearest J 2 neighbor spin interactions. Despite its sim- 
plicity, it is not only in its own right academically inter- 
esting, but also practically provides a window for looking 
inside into the superconductivity mechanism. Recently 
found vanadium-based compounds have provided exper- 
imental realizations of J1-J2 model with J2/J1 ~ 1— 
and again have stimulated research on this charming sys- 
tem. 

As is known, a single layer in undoped cupratc high-T c 
compounds are well described by a square-lattice Heisen- 
berg model displaying long-range antiferromagnetic Neel 
order. The common belief is that upon doping, this long 
range order is destroyed and a different non-magnetic 
phase sets in, which, accompanied by fluctuations, turns 
the system into superconducting phase. One idea sug- 
gests that the effect of doping in destroying the Neel or- 
der might be accounted for by the introduction of frustra- 
tion in original Heisenberg model 9 . Following this sug- 
gestion, many works had been focused on finding such 
phases in frustrated quantum magnet o 10 ' 11 ^ 2 ' 13 . Re- 
cently, the half-filled Hubbard model with both near- 
est neighbor and next-nearest-neighbor hopping term has 
been investigated^. This model in the large U limit is 
equivalent to J1-J2 model. In that paper, variational 
cluster approach shows that d x 2_ y 2-wewe superconduc- 
tivity can also occur at half-filling when the Hubbard 
system is under pressure provided that the frustration 
and the on-site repulsion are not too large. 

However, there is another issue which is worth while 
to be investigated: possible superconducting symmetries 
of J\-J% model under doping. 

At half-filling J\-J 2 model respectively exhibits two 
long-range magnetic orders in two distinct limit: a) the 
Neel phase in the limit of J 2 /J\ — > 0; b) the so called 
collinear phased with ordering wave vector (tt, 0) or 
(0, 7r) in the limit of J2/J1 >> 1. Frustration makes 
the system highly degenerated in the intermediate range 



and may dramatically suppresses the magnetic long range 
orde r 16 ' 17 i 18 ' 19 ' 20 . Hence together with doping, frustra- 
tion may provide the system exotic superconducting sym- 
metries over a broad range of doping and frustration. 

In this work we investigate superconducting symme- 
try of the ground state of doped J\-J 2 model, i.e. t\- 
t-2-Ji-Ji model. To investigate the properties of the 
ground state, a simple yet powerful method is the 
renormalized mean-field theory (RMFT) 21,22 in which 
the kinetic and superexchange energies are renormal- 
ized by different doping-dependent factors gt and g s , 
respectively. Despite of the simplicity of this method, 
it can lead to semi-quantitative even quantitative ex- 
planation of some ground state properties of cuprate 
superconductor a 21 ! 22 ' 23 . In this letter, with the help of 
RMFT, we show that the superconducting symmetry of 
the ground state varies among different types when tun- 
ing the frustration amplitude and doping concentration 
of the system. 

Model — The Hamiltonian of the t\-t 2 -J\-J 2 model 
takes the form of 

H = P d H t P d + H s , 

{nn)a {nnn)<7 

h s = j x J2 Si • Si + J 2 s * ' s j . (!) 

(tin) (nnn) 

Pd = Jl(l — n^riii) is the Gutzwiller projection 

i 

operator— which removes totally the doubly occupied 
states. t\ and t 2 are nearest-neighbor and next-nearest- 
neighbor hopping amplitude. When they are all posi- 
tive, the Hamiltonian represents hole doping case. The 
electron doping case can be achieved via particle-hole 
transformation, changing the sign of t 2 while keeping 
the sign of t\ unchanged 24 . J\ and J 2 are respectively 
the nearest-neighbor and next-nearest-neighbor antifer- 
romagnetic coupling constants, they raise frustration in 
the system. We use t\ as the energy unit and set 
ti/Ji = 3 for conventional reaso n 25 i 26 , J 2 j J\ = {t 2 /ti) 2 
since the superexchanges have relations of J = 4t 2 /U 
with hopping parameters in the large Hubbard U limit. 
And we take 77 = t 2 /ti as frustration amplitude. 
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Method — Renormalized mean-field theory. In the 
frame of RMFT, to investigate the ground state of the 
above mentioned Hamiltonian, the trial state is sug- 
gested to be a projected BCS state \ty) = Pdl^o), where 

l*o) = I\k( u k + u fe4| c1 lfei)l )- And tne projection op- 
erator is taken into account by a set of renormalized 
factors2L22^a, i.e. we have {ip\H\tp) = (i/) \H'\ipo) = 
(ipo\g t H t + g s H s \?po}. In homogenous case the renormal- 
ized factors^. g t = 26/(1 + 6) and g s = 4/(1 + 5) 2 . 

Minimize the quantity W = (H' — fi ^Z i<7 c\ a Ci a )a 
with respect to and i>k, and introduce two mean- 
field parameters A T = (c\^c\ +Tl - 4i c !+tt)° and & = 
5Z (T (cj cr Ci+ T , (7 )o, where r indicates four different bond di- 
rections sketched in FigdJ we get the coupled gap equa- 
tions as follows 

A r = ^-^COBfcrAg/Bg , (2) 
k 

£ r =- N-^cosk^/Eg , (3) 
fc 

where Aj = J2 T At J t /J\ cosfc T , ^ = - 
E r 6r Jr/J\ cos fcr, % = + |A £ | 2 , iV s is the to- 

tal number of sites. In the above equations £^ = 
(-2g t tiJ2 T U/tiCosk T - p)/{\g a J) and £ = fi + 
N~ 1 {^-)o- These gap equations should be solved simul- 
taneously with 5 = N^ 1 J2k£k~/Ek~i th e resulting A's de- 
termine the symmetry of possible superconductivity. 

Although the RMFT cannot provide us a true picture 
of the system in exact half-filling case, the symmetries of 
gap parameters obtained at that point will not change 
under small doping 21 . Thus, in order to investigate pos- 
sible superconducting symmetries of the system under 
doping, our strategy is following: at first, we solve the 
gap equation at half filling, find out the possible sym- 
metries of the mean-field parameters. These symmetry 
states may degenerate at half-filling. Then we switch on 
the doping, compare the energies of different symmetry 
states, we can find the true ground states for different 
doping level and frustration amplitude. 

At half filling, energy if per site is — | g s JJ2k E^, with 
the ansatz that Ek of ground state energy can be written 
as a summation of square terms of cosine functions, after 
some calculation, we find that several symmetry states 
come out, including: (1) D xy + is^-wave state degen- 
erates with D xy + iD a i,-wa,ve state, both of them have 
£,a,b = 0; (2) D a b + iD xy -wa,ve state degenerates with 
D a b + is xy -~w&ve state, £ XtV = in both states. Here i 
means that the difference between the phase <f> a of A a 
and the phase <p x of A x is \<p a — (f> x \ = 7T /2, we use D T1T2 
and s TlT2 to denote <i-wave symmetry and s-wave sym- 
metry on direction t% and ti respectively. 

When changing the doping level and the frustration 
amplitude, these superconducting symmetries compete, 
and energetically, each occupies a specific region in the 
phase diagram, as shown in FigQ] In this phase diagram 
dashed line are boundary of the first order phase transi- 
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FIG. 1: Left picture is schematic structure of tvti-Jx-J^ 
model lattice. We use a, b to distinguish the two different 
orientation diagonal bond. Right picture is phase diagram 
of this lattice, dashed lines are boundary of first order phase 
transition. 
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FIG. 2: Figure(a) shows mean-field parameters of D xy +iD a b 
state as function of S for rj — y/0£. Figure(b) shows param- 
eter functions for r\ = yL5, the inset picture shows phase of 
pairing parameters, its unit is n. 

tion where mean-field parameters change suddenly. The 
bold lines show second order phase transitions. Left pic- 
ture in Fig[T] is a cartoon sketch of the lattice structure 
that we discuss, a and b denote the two different diagonal 
bond. In the following, we are going to discuss the phase 
diagram for positive and the negative £2/^1 case in more 
detail. 

Positive rj cases. — For 77 < \/0.5, as expected, al- 
though the next nearest neighbor interaction frustrates 
the system, it does not destroy the £>-wave supercon- 
ducting symmetry state with A x = — A y and A a .b = 
for all doping level under investigation. For r\ > \/0.5, 
pairing on diagonal links should be considered seriously. 

RMFT calculation shows that when the frustration 
amplitude falls in the range of V0.5 < 77 < V0.7, 
D xy + iD Q f,-wave state with |A X)3/ | > |A 0] &| is the most 
energetically favored state in the small doping region. 
In order to show the dependence on the doping con- 
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FIG. 3: (color online) For r\ = \/0.9, Fig. (a) shows mean-field 
parameters of D xy + iD a b-wa.ve state, while Fig.(b) shows that 
of s xy + i-Dab-wave, the inset of (b) shows the phases of its 
pairing parameters. 



FIG. 4: (color online) Figure(a) shows mean-field parameters 
amplitude of stable state as function of S for 77 = — \/0.6, 
figure(b) shows parameter function of stable state for n = 
— \/1.2, the inset of (a,b) shows the phase of A. 



centration, as an example, we take 77 = -\/0.6 and plot 
D xy + iD a b-waye parameters in Figj^a). With increas- 
ing 6, A b decreases rapidly. One can see from phase 
diagram FigfT]that after S = 0.06, D xy -wa.ve state is en- 
ergetically more favored, and there is a sudden change of 
A a j, in this transition. 

For about V0.8 < 77 < y/l.l in small doping level stable 
state is still D xy + iD a b-~w&ve state, however, in relatively 
high doping level another mixed superconducting sym- 
metry state of s xy + iD a b-w&ve dominates. Again, in or- 
der to illustrate the dependence on doping level, we take 
r] = V0.9 and plot the parameters for D xy + iD a b sym- 
metry and s xy + iD a i, symmetry in Fig[3]^a) and Fig[3fb), 
respectively. Inset of FigG^b) shows the phases of A T in 
different bond directions, and clearly reflects the symme- 
try property of the s^+z-Dafc-wave state. From this inset, 
we can also see that the phases of A a .b is slightly away 
from 7r/2. By comparing the energies of these two differ- 
ent states, the calculation shows that a first order phase 
transition from D xy + iD a b-wave state to s xy + iD a b- 
wave state occurs when increasing the doping level 8 
across 0.06. When 77 is tinily larger than 1, in very small 
doping region stable state most likely is another type of 
D a (, + iD xy -w&ve with |A aj b| > IA^jJ, it occupies really 
an extremely small area and can not be shown explicitly 
in our phase diagram. 

In the region of \f\Ji < 77 < \/2, stable state is fab- 
wave in small doping and will change to s xy + iD a b- 
wave state when increasing the doping level with the 
corresponding phase transition being second order. In 
FigEfb) we show the case of 77 = VT1), it is clear that 
the parameters vary smoothly from Z3 ao -wave type to 
s xy + iD a b-w&ve type when increasing 6. Inset of Figj2]^b) 
clearly shows that in s xy + i£> a b-wave state phase differ- 
ence of A a i, is not exactly n and varies slowly with doping 
concentration. 



For very large 77, it is reasonable to expect a state with 
d-wave pairing only on diagonal bonds, i.e. the state is 
D a h-wave. 

Negative r\ case. — As in the positive 77 case, when |?7| 
takes small value, roughly smaller than V0.5, supercon- 
ducting symmetry of the system is conventional rf-wave 
on plaquette bonds. 

For — V0.7 < 77 < — V0.5, the D xy -w&ve symmetry 
changes into s a b + iD xy -wewe state when increasing dop- 
ing level with the phase transition being first order. It 
should be emphasized that the phase difference between 
A x and A a in the latter state is less than tt/2 and de- 
creases with increasing doping level. FigQJa) shows that 
for 77 = — V0.6, when <5 > 0.16, D xy -wa,ve state discontin- 
uously changes to s a & + iD xy -wave state. Inset shows the 
phase of pairing parameters, at the critical point symme- 
try of superconductivity changes suddenly. 

When — VL5 < 77 < — vLl, Z) a b-wave is the stable 
state in high doping level while in small doping level sta- 
ble state is s xy + i-D a &-wave symmetry. 77 = — y/1,2 case 
is shown explicitly in Figj4jb) as an example. With in- 
creasing S pairing on diagonal bonds A x ^ y decrease to 
zero rapidly and s xy + i-D ao -wave state varies to £> a fc- 
wave smoothly through a second order phase transition. 
-Dab-wave state occupies more range of 6 when increasing 
1 77 1 till 1 77 1 = \/F5, after that D a b is the only stable one 
for all doping level. 

For — vLl < 77 < — V0.7, in small doping level 
D xy + iD a b state has more energy gain while s xy + iD a b 
is more stable at higher doping level. FigOJa) shows how 
parameter amplitudes varied with doping concentration 
at 77 = — Vol) for best energy stable state and FigJSfb) 
shows the phase changing of stable state. The sudden 
changes of the parameter amplitudes and phases clearly 
indicate that the phase transition here is first order. 

Summary and discussion — In summary, we have inves- 
tigated the possible superconducting symmetries of frus- 
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FIG. 5: (color online) Figure(a) shows mean-field parameters 
amplitude of stable state as function of S for n = — \/0.9, 
figure(b) shows phase of pairing parameters. At the critical 
point symmetry of superconductivity change for D xy + %D a ^ 
to s xy + iD a b wave state. 



trated t\-t2-J\-Ji model by RMFT method. In terms 
of frustration amplitude t%jt\ and doping concentration 
S a phase diagram for possible superconducting symme- 
try states has been presented. We have shown that in 
weakly and strongly frustrated cases pure d-wave states 
are on plaquette or diagonal bonds, respectively. When 
t\ and t2 arc comparable, pairing parameters on plaque- 
tte bonds and diagonal bonds are both finite and mixed 
symmetry states such as D + ID, D + is appear. Each 
symmetry state occupies a specific region in the phase di- 
agram, transitions between these states are the first order 
or second order, the latter one corresponding to relatively 
large J 2- The two kinds of time-reversal broken mixed 
states are consistent with those found in two-dimensional 
Fermi liquid with attractive interaction in both s and D 
channel 31 , D + is time- reversal broken state may also ap- 
pear in the pseudogap phase 32 . We hope that our calcula- 
tion may shed some new light on the study of mechanism 
of superconducting order in high-T c . 
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